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Abstract 

A popular three-dimensional reduction of the Bethe-Salpeter formalism for the description 
of bound states in quantum field theory is the Salpeter equation, derived by assuming both 
instantaneous interactions and free propagation of all bound-state constituents. Numerical 
(variational) studies of the Salpeter equation with confining interaction, however, observed 
specific instabilities of the solutions, likely related to the Klein paradox and rendering (part 
of the) bound states unstable. An analytic investigation of the problem by a comprehensive 
spectral analysis is feasible for the reduced Salpeter equation with only harmonic-oscillator 
confining interactions. There we are able to prove rigorously that the bound-state solutions 
correspond to real discrete spectra bounded from below and are thus free of all instabilities. 
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1 Introduction 



The most widely explored three-dimensional reduction of the Bethe-Salpeter formalism [1] 
for the description of bound states within quantum field theory is the Salpeter equation [2]. 
The Salpeter equation controls the Salpeter amplitude, which in momentum space encodes 
the distribution of relative momenta of all bound-state constituents. In elementary particle 
physics its application to quantum electrodynamics (QED) and quantum chromodynamics 
(QCD) has met considerable success. In particular, within the latter realm it has evolved to 
a well-established standard tool for describing from first principles hadrons as bound states 
of quarks, confined by the strong interactions. Surprisingly or not, however, the solutions of 
Salpeter's equation with confining interactions have been numerically shown to develop for 
some Lorentz structures of the Bethe-Salpeter kernel representing all interactions between 
the bound-state constituents instabilities which cause states expected to be stable to decay. 

In contrast, the reduced Salpeter equation [3-7], derived from the full Salpeter equation 
by neglecting some of the interaction terms (Sec. 2), offers the chance to study the question 
of stability analytically. Such instabilities should arise first for pseudoscalar states (Sec. 3); 
there, stripping off all angular variables simplifies [8-15] the reduced Salpeter equation to a 
single integral equation (Sec. 4). Harmonic-oscillator confining interactions (Sec. 5) have a 
big advantage: In momentum space, they are represented by a simple Laplacian, converting 
thus our integral to differential equations (Sec. 6). Their analysis proves that, for all famous 
Lorentz structures, including one proposed by Bohm, Joos, and Krammer (B JK hereafter) 
[16], studied in Ref. [17] and more recently used, among others, by a group in Bonn [18-21], 
any bound-state solution is related to a real (Sec. 7) and discrete (Sec. 8) energy eigenvalue, 
and thus stable. Similar considerations can be applied to the full Salpeter equation (Sec. 9). 



2 Bethe— Salpeter Formalism in Instantaneous Limit 



We are interested in instantaneous approximations to Bethe-Salpeter equations describing 
bound states composed of a fermion and an antifermion. Let these constituents, denoted by 
% — 1,2, carry the momenta p\ and P2, which will enter also in terms of the total momentum 
P = P1+P2 and, for a real parameter 77 or (, the relative momentum p = p\—r\ P = ( P—p2- 

2.1 Salpeter equation 

Any derivation of the Salpeter equation [2] as one of a variety of possible three-dimensional 
reductions of the Bethe-Salpeter formalism is based on just two fundamental assumptions: 
First, one obtains the instantaneous Bethe-Salpeter equation if all interactions between the 
bound-state constituents are instantaneous in the center-of-momentum frame of the bound 
state. In the Bethe-Salpeter equation [1] all interactions between bound-state constituents 
are encoded in their integral kernel, K. The instantaneous approximation then implies that 
K depends only on the spatial components of the relative momenta involved: K = K(p, q). 
Second, one arrives at the Salpeter equation if every bound-state constituent propagates as 
free particle with effective mass, m. Salpeter's equation adopts the free fermion propagator 



S (p,m) 



p' + m 



e I . 



p 1 — m + is 



V 



> 2 — m 2 + i e 
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Let us present those aspects of the Salpeter equation which will be relevant for stability. 
We introduce the one-particle energy E^p), the one-particle Dirac Hamiltonian H^p), and 
the energy projection operators Af (p) for positive or negative energy of particle % — 1, 2 by 



Ei(p) = ^P 2 + m*, 2 = 1,2, 
Hiip) = -foil-P + rrii), i = 1,2 , 

related, for i = 1, 2, by Hf(p) = Ef(p) and Af (p) H t (p) = H^p) Af (p) = ± E^p) Af (p), 
all operators Af (p) fulfilling Af (p) Af (p) = Af (p), Af (p) Af (p) = 0, Af (p)+Ar(p) = 1. 
In terms of the above abbreviations, the Salpeter equation [2] , for bound states of a fermion 
(of mass mi and momentum p-J and an antifermion (of mass m 2 and momentum p 2 ), reads 

<$>( n \ f ( Af (Pi) 7o [K{p, q) $(g)] 7o A 2 (p 2 ) 

W " 7 (2tt)3 v P - El ( Pl ) - E 2 (p 2 ) 

Mgi) 7o g(g, g) 7o Af (p 2 ) \ 
P + ^i(Pi)+^(p 2 ) /' 
Accordingly, every solution $(p) of the Salpeter equation has to satisfy the two constraints 

Af $(p) A+(p 2 ) = Ar( Pl ) $(p) A 2 -(p 2 ) = , 

which, by considering their sum or difference, prove to be equivalent to a single constraint: 

or, equivalently, 

With the help of the decomposition of unity in terms of projection operators Af (p), that is, 
1(8)1 = Af (p 1 )®A+(p 2 )+A+(p 1 )®A 2 ^(p 2 )+Ar(p 1 )®A+(p 2 )+Ar(p 1 )®A 2 (p 2 ), a Salpeter 
amplitude $(p) may be, in general, decomposed into the components Af(p 1 ) $(p) Af (p 2 ): 

$(p) = Af(p 1 )$(p)A+(p 2 )+Af(p 1 )$(p)A 2 (p 2 ) 
+ Ar(Px) $(p) A+(p 2 ) + Ar(Px) $(p) A 2 (p 2 ) . 

The above constraint (s) on the solutions $(p) of the Salpeter equation (1), arising from its 
specific projector structure, halve, in fact, the number of independent components of $(p): 

$(p) = Af $(p) A 2 (p 2 ) + Ar( Pl ) $(p) A+(p 2 ) . (3) 

Assuming the Lorentz structures of the effective couplings of the bound fermions % = 1, 2 to 
the corresponding interaction potentials to be identical, the Bethe-Salpeter kernel K(p, q) 
is, quite generally, the sum of products of a tensor product Y®Y of generic Dirac matrices Y 
and a Lorentz-scalar associated interaction function Vr(p, q)' K(p, q) = J2r Vr(p, q) Y®Y. 
More precisely, the action of the kernel K(p, q) on the Salpeter amplitude $(p) is given by 

[K(p,q)^(q)]=Y / Vr(p,q)Y^q)Y . 
r 
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2.2 Reduced Salpeter equation 

The projection operators Af(p), i = 1, 2, for positive or negative energy satisfy the identity 



Af(p) 



C7" 1 Af(p)C7 =A?(p), 



1,2, 



where C denotes the usual Dirac-space charge-conjugation matrix. Accordingly, the second 
term on the right-hand side of the Salpeter equation, Eq. (1), is that interaction term which 
is related to the negative- energy components Aj" (p-J $(p) A^(j> 2 ) = Af(Pi) 3>(p) [Aj (p 2 )] c 
of the Salpeter amplitude $(p). Assuming that this term's contribution may be reasonably 
ignored relative to that of the first term yields the so-called reduced Salpeter equation [3-7] 



[P -E 1 (p 1 )-E 2 (p 2 )]$(p) 



d 3 q 



At( Pl ) lo [K(p,q)<S>(q)] lo A 2 (p 2 ) . (4) 



This neglect might be justifiable for nonrelativistic and weakly bound systems composed of 
heavy constituents, such that, on the average, P — E 1 (p 1 )— E 2 (p 2 ) <C Po+E^p^+E^p^. 
More rigorously, one would like to be sure, at least, that, for appropriate expectation values, 

1 1 



P + E 1 (p 1 ) + E 2 (p 2 ) P - E 1 (p l ) - E 2 (p 2 ) ■ 

Formally, this reduction of the (full) Salpeter equation (1) to the reduced Salpeter equation 
(4) can be accomplished by subjecting the Salpeter amplitude $(p) to any of the equalities 



A-( Pl )<S>(p) = , $(p)A+(p 2 )=0, 



or, equivalently, 

ififo) $(p) = E 1 (p 1 ) $(p) , $(p) H 2 (p 2 ) = -E 2 (p 2 ) $(p) ; 

because of the particular projector structure of the resulting reduced Salpeter equation (4), 
imposition of one of these two constraints automatically implies the other for the solutions. 
These equalities entail the constraints Af(p 1 ) $(p) A2~(p 2 ) = Aj"(p 1 ) $(p) A2"(p 2 ) = also 
satisfied by each solution of the Salpeter equation (1), as well as, in addition, the constraint 
Aj"(p 1 ) $(p) A^ (p 2 ) = 0. As trivial consequence, any solution $(p) of the reduced Salpeter 
equation (4) is necessarily of the unique component structure $(p) = Af(p 1 ) $(p) A2"(p 2 ). 



3 Pseudoscalar Bound States 

The states easiest to investigate are bound states composed of a fermion and the associated 
antifermion, which guarantees a well-defined behaviour with respect to charge conjugation. 
The masses of a particle and the corresponding antiparticle are, of course, identical, that is, 
m i — m 2 — m - Now, in the center- of -momentum (or rest) frame of the two-particle system 
under study, defined by P = 0, which implies p = p 1 = — p 2 , the time component Pq of the 
total momentum P reduces to the bound-state mass eigenvalue M, i.e., Po = M. All indices 
i — 1,2, distinguishing the two bound-state constituents, may then be dropped throughout 
the following analysis: E ± (p) = E 2 {p) = E{p) = E(p) = ^Jp 2 + m 2 , where p = |p| = \fp z ] 
and, mutatis mutandis, for the Hamiltonians Hi(p) and energy projection operators Af (p). 
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For the sake of definiteness we focus, in what follows, to the case of fermion-antifermion 
bound states of total spin J, parity P = (— 1) J+1 , and charge-conjugation quantum number 
C = (— 1) J (which entails C P = — 1 for all J). In usual spectroscopic notation, these states 
are denoted by n 1 Jj (n = 1, 2, 3, . . .). As the perhaps simplest example within this context, 
precisely these physical systems have been studied rather frequently (cf., e.g., Refs. [8-15]). 

More specifically, we consider, for reasons of simplicity, bound states of spin J = 0, that 
is, pseudoscalar bound states, with spin-parity-charge conjugation assignment J PC = CT + , 
spectroscopically called 1 So- In nature such systems are realized and observed, for instance, 
in the realm of quark-ant iquark bound states, in form of the pion and its radial excitations. 

In the case of the Salpeter equation (1), as a consequence of the constraints discussed in 
Subsec. 2.1 the most general expansion of the Salpeter amplitude $(p), over a complete set 
of Dirac matrices, involves not the full 16 but only eight independent Salpeter components. 
For the description of 1 Jj states, only two of the latter, called 0i(p) and 2 (p), are relevant. 
With our notation for one-particle energy E(p) and Dirac Hamiltonian H(p) introduced in 
Subsec. 2.1, full 1 Jj Salpeter amplitudes $(p) for fermion and antifermion of equal mass m 
and internal momentum p thus read, in the center-of-momentum frame of the bound state, 



Hp) 



75 



In the case of the reduced Salpeter equation (4), the one additional, two-faced constraint 
analyzed in Subsec. 2.2 entails, possibly taking into account 75 H(—p) = —H(p) 75, for the 
two independent Salpeter components 0i(p) and 2 (p) in the 1 Jj Salpeter amplitude $(p), 
0i (p) = 02 (p) = 0(p)- Consequently, the generic 1 Jj reduced-Salpeter solutions <3>(p) read 

Hp) = 0(p) H{p) ^ f {p) 75 = 2 0(p) A+b) 75 • 

E(p) 



4 Radial Eigenvalue Equations 

For any (instantaneous) Bethe-Salpeter interaction kernel K(p, q) of convolution type, i.e., 
K{p, q) — K(p—<l) an d therefore Vr(p, q) = Vr(p—q), by factorizing off all dependence on 
angular variables encoded in corresponding (vector) spherical harmonics both the Salpeter 
equation (1) [8,10] and its reduced version (4) [11] can be converted into equivalent systems 
of coupled equations for the radial factors of all relevant independent Salpeter components. 
For a fixed Lorentz structure of the kernel, the interactions experienced by the bound-state 
constituents enter in such set of equations in the form of Fourier-Bessel transforms Vl{p, q) 
(L — 0, 1, 2, . . .) of some spherically symmetric static potential V(r) in configuration space: 

00 

VL(p,q) = 8nJdrr 2 j L (pr)j L (qr)V(r), L = 0,1,2,... , 


where j n (z), for n = 0, ±1, ±2, . . . , label the spherical Bessel functions of the first kind [22]. 

According to Sec. 3, for pseudoscalar fermion-antifermion bound states each solution of 
the reduced Salpeter equation (4) involves only one independent Salpeter component, 4>{p). 
Consequently, the aforementioned system of radial equations collapses to a single equation. 
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The application of the radial reduction to the reduced Salpeter equation (4) then yields the 
radial eigenvalue equations, for interactions of Lorentz-scalar Dirac structure, r<g>r = 1(8)1, 



2E(p)<j>(p) 



1 f dqq 2 



2 J (2nf 



1 + 



m 



E{p)E(q)_ 



Vo(p,q) - 



pq 



E(p) E(q) 



Vi(p,q) 



4>{q) = M(f ) {p) 



for interactions of time-component Lorentz- vector Dirac structure, r <g) T — 7 <g) 7 , 



2E(p)</>(p) + -J 



1 f dqq 2 



{2nf 



1 + 



m 



E(p)E(q), 



Vo(p,q) + 



pq 



E(p) E(q) 



Vi(p,q) 



<P(q) = M<P(p) , 



for interactions of Lorentz-vector Dirac structure, r ® V = 7^ ® 7^, 



2E(p) ( f>{ P ) 



2 - 



m 



dqq 2 

Wj*V E(p)E(q) : 



Vo(p,q)<j>(q)=M<j>(p) , 



for interactions of Lorentz-pseudoscalar Dirac structure, r <8> T = 75 ® 75, 



2E{p)(t>{p)- 



1 [dqq 7 



2 J (2tt)' 



1 - 



Vb(p,9) - 



K(P,?) 



0(g) = M0(p) 



E(p)E(q)J uv "'^ £(p)£(g) 

and, for interactions of BJK [16,17] Dirac structure, r®r = | (7^ <8> 7 M + 75 (g> 75 — 1 <8> 1) , 

v1 j 2 
dqq 



2E(p)<P(p) + 



V (p,q)(f>(q) = M ( f>(p) 



For any mass eigenvalue M such that M—2 E(p) 7^ the above radial eigenvalue equations 
are all of the form of a homogeneous linear (Fredholm) integral equation of the second kind. 



5 Harmonic- Oscillator Interaction 

The reduced Salpeter equations of Sec. 4 may be discussed to a large extent analytically by 
focusing to harmonic-oscillator interactions described by the configuration-space potential 

V(r) = ar 2 , a = a* 7^ , r = \x\ . 

For this choice of V(r) we are able to determine analytically all potential functions Vl(p, q) 
entering in the radial eigenvalue equations, by taking advantage of the differential equation 
[22] satisfied by all spherical Bessel functions, generically called w n (z) (n = 0, ±1, ±2, . . .): 

d 2 d 
z 2 —w n (z) + 2z —w n (z) + [z 2 - n (n + 1)] w n (z) = . 
dz A dz 

Defining, as radial relic of the Laplacian A = V- V, the second-order differential operators 

J + 2 d _ LJL + V L = QA ^_ 
v dj9 z p dp p A 
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the spherical Bessel functions of the first kind, jL{p r ), in the potentials Vl{p, q) thus satisfy 

D ( p L) J L (pr) = -r 2 j L (pr) , L = 0,1,2,... . 

This relation allows to replace the harmonic-oscillator potential r 2 in the potential function 
q) by the differential operator D^ L \ Hence, by means of the "orthogonality relations" 



drr 2 j L (pr) j L (qr) 



7T 

2/7 



6(p-q), L = 0,1,2, 



involving Dirac's delta distribution the potential functions for harmonic oscillators become 

(2ix) 2 a 



v L {p,q) = 



D^S(p-q) , L = 0,1,2,... . 



(5) 



6 Ordinary Differential Equations (of Second Order) 

For the potential functions (5) representing some harmonic-oscillator interaction, all radial 
integral eigenvalue equations derived in Sec. 4 simplify to second-order homogeneous linear 
differential equations; these read, for kernels of Lorentz-scalar Dirac structure T®T = 1<8>1, 



2E(p) + a 



E\p) 



^m 2 (p 2 — 5m 2 ) 2m 2 p d 



or 



2E{p) + a 



2E\p) 

'2p 2 + 3 m 2 
, lE\p) 



D (o) 

E\p) dp E*(p) p , 



4>(p) = M 4>(p) 



+ 



m 



E(p) p E(p) : 



<j>(p) = M0(p) 



for kernels of time- component Lorentz-vector Dirac structure r (8 T = 7 (8 7° 



2E(p) + a 



'2p 2 + 3m 2 



2 E\p) 

for kernels of Lorentz-vector Dirac structure T 



4° } 



4>{p) = M 4>{p) , 



(6) 



(7) 



3am 4 2am 2 p d / 



m 



2 



i- 1 



E 2 {p) J p 



4>{p) = M<t>(jp) 



or 



2 E(p) + a 



( JH— £>(o) _i 2 D (0) " 



(f>(p) — M<f)(p) , 



for kernels of Lorentz-pseudoscalar Dirac structure T (8 T = 75 (8 75, 



2E(p) + a 



2p 2 + 3m 2 
2^ 4 (p) 



<j>(jp) = M <j>{p) , 



(8) 



(9) 



and, for kernels of the BJK [16,17] Dirac structure r<8r = \ (7^ (8 7 M + 75 <E> 75 — 1 <E> 1) , 



2£(p) - aDf ] \ <p(p) = M<p(p) . 



(10) 



Note that, for the Lorentz pseudoscalar r<8T = 75(875, the reduced Salpeter equation with 
harmonic-oscillator potential is represented by a pure multiplication operator. This implies 
that the resulting spectrum is purely continuous. That is, there are no bound states at all. 
This fact is presumably not evident from the general representation of the reduced Salpeter 
equation with pseudoscalar Lorentz structure as radial integral equation, as given in Sec. 4. 
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6.1 Transformation to a (zero-eigenvalue) Schrodinger equation 

For m^fl, the ordinary differential equations that represent the reduced Salpeter equation 
(4) for a harmonic-oscillator interaction of Lorentz- scalar or Lorentz-vector Dirac structure 
do not (yet) resemble the familiar form of Schrodinger eigenvalue equations. However, both 
of them can be easily cast into the form of the ordinary differential equation of second order 



d 2 , * d , . . 

.--2 g ( P )- + h( P ) 



[ID 



involving two given functions, g(p) and h(p) . Here, similarly to the role of the mass m of the 
two bound-state constituents and the coupling a of the harmonic-oscillator interaction, the 
bound-state mass eigenvalue, M, enters — in the function h(p) only — as a parameter. Then, 
by substitution of the amplitudes (f>(p) by 4>{p) = f(p) ip(p), with the transforming function 



f(p) =pexp 



dpg(p) 



d 
dp 



g(p) 



fip) 



p 



determined up to an irrelevant constant, Eq. (11) becomes the eigenvalue equation for ipip) 

d 2 2d 



+ U(p) 



4>(p) = , 



dp 2 p dp 

corresponding to eigenvalue 0, of the Schrodinger operator — D^+U (p), with the potential 

1 d 2 f 



U(p) = h(p) 



9 (P) 



;i2) 



fip) d P 2 

here the second equality follows most easily from the differential equation satisfied by fip). 
Let us apply this transformation to the two cumbersome Lorentz structures under concern. 

• In the case of a Lorentz- scalar Bethe-Salpeter kernel, Y®Y = 1<8>1, the integration of 



g{p) = - 
p 



p 



E 2 (p) 



gives f(p) = E{p) [as may be guessed from Eq. (6)] whereas the potential U{p) reads 

mw-iEm-^—L 



U{p) 



am* 



m* 



(13) 



2E 2 (p) ' 

• For the case of a Lorentz-vector interaction kernel, T®T = 7^(87^, the integration of 



1 



m 2 p 



9 ^ p + E 2 (p) [E 2 (p) 



p z 



(which can be performed most conveniently by partial fraction decomposition) yields 

ftp) = E & ■ 
^E 2 {p)+p 2 ' 

reading off h(p) and inserting it into the expression (12) for the potential U(p) entails 

E 2 (p) [2 E(p) — M] 2m 2 p 2 



U(p) = 



a[E 2 (p)+p 2 ] E 2 (p)[E 2 (p)+p 2 } 2 • 



(14) 



With these formulations of harmonic-oscillator reduced Salpeter equations at our disposal, 
we may take advantage, in Sec. 8, of some general results derived for Schrodinger operators. 
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7 Self-Adjointness of "Reduced-Salpeter" Operators 



The main goal of our investigation is the thorough analysis of the qualitative features of the 
spectra of the reduced Salpeter equation (4) with harmonic-oscillator interaction of various 
Lorentz structures by close inspection of the resulting differential equations given in Sec. 6. 

In this respect, the first problem that arises is the question of the self-adjointness of the 
operators in our (differential) equations (6) through (10). It is straightforward to show that 
all these operators are self-adjoint: multiplication by a real-valued function clearly defines a 
self-adjoint operator and by integration by parts it is easy to convince oneself that and 



m 



1 £K°) 
E\p) * 



2p d\_ m 2 (0) 1 
E\p)dp) E{p) p E(p) + 



3 m 4 

EHp) 



or 



m 



£>(0)_L_ 

E{p) p E(p) 



m 2 2m 2 p d 3 m 4 

E 2 {p) p ~ E 4 (p) dp~ ~ E e (p) 



also represent self-adjoint operators. This implies that the corresponding spectra are real. 

The reality of all eigenvalues M can be also inferred from Eqs. (18) and (16) of Ref. [11] 
along the lines of argument presented (for the case of the full Salpeter equation, however) 
in Sec. II of Ref. [8] or Sec. II of Ref. [23]. From Eq. (4), all Salpeter amplitudes $(p) satisfy 



M 



1 



d 3 p 
J2mf 



Tr 



*t(p)$(p) 



+ 



d 3 p 

d 3 p 
(2tt) 3 J (2tt 



[E 1 (p) + E 2 (p)]Tr[^(p)^(p) 
d 3 q 



^\/r(p,g)Tr[$t(p) 7o r$( q )r7o 



In this relation, both the integral multiplied by M on the left-hand side, that is, the "norm" 



l$l 



d 3 p 



Tr 



&(p)$(p) 



of each Salpeter amplitude $(p) emerging from the reduced Salpeter equation, and the first 
term on the right-hand side are certainly real; the second term on the right-hand side is real 
provided all Lorentz-scalar potential functions Vr(p, q) satisfy V^(q, p) = Vr(p, q) and the 
Dirac matrices Y satisfy 70 1^ 70 = ±I\ which implies that the matrices r = 7o V are (anti-) 
Hermitian, i.e., H = ±f. In contrast to the full Salpeter equation, for the reduced Salpeter 
equation the norm ||$|| 2 of all nonzero solutions $(p) is definitely nonvanishing. Therefore, 
all mass eigenvalues M are guaranteed to be real for reasonable interaction kernels K(p, q). 



8 Spectra and Stability of the Bound States 
8.1 General considerations 

Now let us investigate in turn the spectra corresponding to the different Lorentz structures. 
Our task is considerably facilitated by making use of a fundamental theorem [24] about the 
spectra of Hamiltonians with potentials increasing without bounds: a Schrodinger operator 
H = — A+ V, defined as a sum of quadratic forms, with locally bounded, positive, infinitely 
rising potential V(x), that is, V(x) — > 00 for \x\ — > 00, may be shown [24], by application of 
the well-known minimum-maximum principle [24-26] , to have a purely discrete spectrum. 
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In the case of harmonic-oscillator interactions of Lorentz- scalar or Lorentz-vector Dirac 
structure, we did not succeed, for m ^ 0, to reformulate the radial eigenvalue equations (6) 
and (8), which fix the bound-state mass eigenvalues M, as standard Schrodinger equations. 
However, in Subsec. 6.1 we managed to cast such unruly eigenvalue equations into the form 

[-D<® + U(p)]i/,(p)=0, 

with U(p) given in terms of well-defined potentials Ui(p), U 2 (p) by U (p) = U 1 (p)+M U 2 (p). 

Accordingly, let us analyze the Hamiltonian operator Hu = —A+U. For U = Ui+M U 2 
satisfying the assumptions of the above "infinitely-rising-potential theorem," the spectrum 
of this operator Hu is, for any M, entirely discrete. In other words, it consists exclusively of 
isolated eigenvalues e(M) of finite multiplicity that depend, of course, on the parameter M. 
Every single zero of these functions £j(M) (i G Z) defines a bound-state mass eigenvalue M 
of Eq. (6) or Eq. (8). The derivative of each such function £j(M) with respect to M is given, 
in accordance with the Hellmann-Feynman theorem [27], by the expectation value over the 
associated eigenstate \i) = 1) of the derivative of this operator Hu with respect to M: 



dM 



SHr 



(i\U 2 \i) 



dM 

Then, if for all eigenvalues £j(M), this derivative is strictly definite, i.e., if for given i either 

dei 



dM 

or 



(M) > V M 



^(M)<0 VI 
dM y 1 

holds, the discreteness of the spectrum of the Hamiltonian Hu for appropriate potentials U 
translates into the discreteness of all eigenvalues M of the bound-state equations (6) or (8). 

8.2 Lorentz-scalar kernel: r <g> T = 1 <g> 1 

For massless bound-state constituents, that is, for m — 0, our harmonic-oscillator reduced 
Salpeter equation with Lorentz-scalar interaction kernel Y®Y — 1(8)1, Eq. (6), simplifies to 

= Mcf)(p) . 

This relation involves a pure multiplication operator, which possesses a purely continuous 
spectrum but no eigenvalues at all. Consequently, the harmonic-oscillator reduced Salpeter 
equation with Lorentz-scalar kernel will not describe bound states of massless constituents. 

In the case of nonvanishing masses m ^ of the bound-state constituents, the auxiliary 
potential U{p) of Eq. (13) satisfies, for negative harmonic-oscillator coupling a = — \a\ < 0, 
all requirements of the "infinitely-rising-potential theorem." The term U 2 multiplied by M, 

E\p) 




U 2 (p) 



,2 



am* 

is, for a < 0, obviously negative definite. Hence, all our considerations of Subsec. 8.1 apply: 

dM a 

Thus, for nonzero bound-state constituents' mass the resulting spectrum is purely discrete. 
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Figure 1: First, lowest-lying eigenvalues £j(M), i — 0, 1, . . . , 9, of the auxiliary Hamiltonian 
Hjj = —A+U with effective potential U of Eq. (13), corresponding to the reduced Salpeter 
equation with harmonic-oscillator interaction of Lorentz-scalar structure, T®r = 1(8)1, for 
bound-state constituents' mass m = 1 and a "binding" coupling a = —10 (arbitrary units). 



Fig. 1 depicts, for a Lorentz-scalar kernel, the typical behaviour of the discrete auxiliary 
eigenvalues e(M) for negative harmonic-oscillator couplings a < producing bound states. 



8.3 Time- component Lorentz-vector kernel: r <S> T = 7 (g> 7 

For arbitrary mass m > of the bound-state constituents, our harmonic-oscillator reduced 
Salpeter equation with a time-component Lorentz-vector interaction kernel TcgF = 7°£g>7°, 
Eq. (7), yields a Schrodinger equation with an effective potential V(p) in momentum space. 
For massless particles, i.e., for m — 0, this potential V(p) is singular at p = (Fig. 2): 

vm = ^ + V 

a p z 

This potential, however, is the sum V{p) = W(p) + l/p 2 of a linear potential W{p) = 2p/a, 
which is positive for a positive slope of this linear rise, that is, for all a > 0, and the singular 
but positive function 1/p 2 . Thus, the functions V(p) and W(p) are related by the inequality 
V(p) > W(p). A suitable combination [28-32] of the minimum-maximum principle [24-26] 
with the resulting operator inequality Hy = — A+V^ > H w = — A+W allows to show that 
any discrete eigenvalue of Hy is bounded from below by a corresponding eigenvalue of H w . 
Since, by the "infinitely-rising-potential theorem," the spectrum of Hw is purely discrete- 
see our detailed discussion in Subsec. 8.4 — , Hy necessarily has a purely discrete spectrum. 
Of course, the same result is obtained by suitable generalization of the theorem of Ref. [24]. 
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(b) 

Figure 2: Qualitative behaviour of the effective potential, V(p), in the differential-equation 
representation (7) of the reduced Salpeter equation with harmonic-oscillator interaction in 
time-component Lorentz-vector kernel, r®r = 7°<S>7°, for a = 10 and (a) vanishing mass 
m = and (b) nonvanishing mass m = 1 of the bound-state constituents (arbitrary units). 
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For nonvanishing mass of the bound-state constituents, m/0, the effective potential is 
fully compatible with the assumptions of the "infinitely-rising-potential theorem" (Fig. 2): 

■ . 2E(p) 2p 2 + 3m 2 . 

v ( p »=^ + ^r (a>0) - 

Here, too, our inevitable conclusion is that the emerging spectrum must be purely discrete. 



8.4 Lorentz-vector kernel: r (g) Y = 7^ (g) 7^ 

For massless bound-state constituents, that is, for m — 0, our harmonic-oscillator reduced 
Salpeter equation with Lorentz-vector interaction kernel r<g>r = 7 M £g>7 M , Eq. (8), becomes 

2 (p - a D®) <j>{p) = M<f)(p) . 

Of course, this is nothing but the standard nonrelativistic Schrodinger equation with linear 
potential, V{p) = p/a, which is equivalent to the Airy differential equation [22]. To see this, 
introduce a reduced radial wave function ip{p) = p<p(p) and perform the change of variables 

1 / M\ 

to a convenient dimensionless variable, z, in order to arrive at the Airy differential equation 

■j-w(z) = zw(z) 

for w(z) = (p(p). Thus, for m = and a positive slope of this linear potential, i.e., for a > 0, 
the "infinitely-rising-potential theorem" guarantees the pure discreteness of this spectrum. 

For nonvanishing bound-state constituents' masses, i.e., m ^ 0, we again have to invoke 
the transformation to an auxiliary Hamiltonian Hy performed in Subsec. 6.1. The resulting 
effective potential U (p) of Eq. (14) is, for positive harmonic-oscillator couplings a > 0, fully 
compatible with all the needs of the "infinitely-rising-potential theorem." Its "M-part" U 2 , 

E 2 {p) 

U ^ P) = ~a[E 2 (p)+p 2 } ' 

is, for all a > 0, negative definite: U 2 < 0. Following our line of argument of Subsec. 8.1, one 
finds that, also for nonzero bound-state constituents' mass, the spectrum is purely discrete. 
Without surprise, the behaviour of the discrete auxiliary eigenvalues e(M) for a > is for a 
Lorentz-vector kernel rather similar to that of their Lorentz-scalar counterparts (cf. Fig. 3). 



8.5 Lorentz-pseudoscalar kernel: r (g) T = 75 (g) 75 

Already in Sec. 6, our inspection of the differential equations revealed that for any kernel of 
Lorentz-pseudoscalar type the spectrum is purely continuous: there is no stability problem. 



8.6 BJK Lorentz structure: r <g> V = | (7^ <g> 7^ + 75 <g> 75 - 1 <g> 1) 

For all m > 0, a harmonic-oscillator reduced Salpeter equation (10) for a kernel of the BJK 
Lorentz structure r<g>r = \ (7^ ® 7 M + 75 ® 75 — 1 ® 1) is just the eigenvalue equation of a 
Schrodinger operator with potential V(p) = 2E(p)/a. From the "infinitely-rising-potential 
theorem" we safely conclude that the spectrum of this operator is purely discrete for a > 0. 
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Figure 3: First, lowest-lying eigenvalues £j(M), i — 0, 1, . . . , 9, of the auxiliary Hamiltonian 
Hu = —A+U with effective potential U of Eq. (14), corresponding to the reduced Salpeter 
equation with a harmonic-oscillator interaction of Lorentz- vector structure T^)T = 7^£g>7 M , 
for bound-state constituents' mass m = 1 and "binding" coupling a = 10 (arbitrary units). 



9 Digression: Investigation of Full Salpeter Equation 

It goes without saying that analogous spectral analyses [33,34] can be envisaged [35] for the 
Salpeter equation (1). There, however, we expect to face more severe problems because any 
solution $(p) of the Salpeter equation (1) involves more than one independent component. 
Thus a Salpeter equation with harmonic-oscillator interaction yields a system of more than 
one second-order differential equation or, equivalently, a higher-order differential equation. 

Hence, let us recall some general features of the eigenvalues of the full Salpeter equation 
(1), regarded as an eigenvalue equation for the Salpeter amplitude $(p); these observations 
emerge either from the relationship, that is, more precisely, the equivalence, of the Salpeter 
equation to the well-known random phase approximation (RPA) familiar from the study of 
collective excitations in nuclear physics or from the inspection of the Salpeter equation (1). 
The normalization condition for Bethe-Salpeter amplitudes arising as a consequence of the 
inhomogeneous Bethe-Salpeter equation [1] motivates the introduction of a "norm" ||$|| of 
the Salpeter amplitude [8,10,23,36], given by (cf. Eq. (2.9) of Ref. [8] or Eq. (9) of Ref. [10]) 



l$l 



d 3 p 
(2^)3 



Tr 



$t(p) 



' ffi(p) 

Mp) 



Hp) - Hp] 



H2{-P) 

Mp) 



Tr 



•»<P)f{g*<P) 



14 



The two expressions on the right-hand side of the above definition are, of course, equivalent 
by virtue of the constraint (2) satisfied by any solution of the Salpeter equation. Because of 
the Hermiticity Hj(p) = Hi(p) of the single-particle Dirac Hamiltonian Hi(p),i = 1, 2, the 
square ||$|| 2 of this "norm" is certainly real. It is, however, not necessarily positive definite. 
(Hence, the Salpeter "norm" ||<3>|| is not necessarily real, whence the quotation marks.) The 
spectrum of Salpeter 's equation then exhibits the following characteristics [8,10,23,37-40]. 

• The Salpeter equation can be shown to be of the same algebraic structure as the RPA 
equation [23,37-40]. In other words, Salpeter's equation and the RPA equation prove 
to be equivalent [38-40] . Now, any RPA equation can always be rewritten in the form 
of a self-adjoint eigenvalue equation for the square of the energy [39,40]. Accordingly, 
the square of the energy is guaranteed to be real, which implies that all eigenvalues of 
such RPA equation are either real or purely imaginary [37]. For the Salpeter equation 
(1) this means that the squares Pq of the energy eigenvalues Pq of a given bound state 
are real or (in the center-of-momentum or "rest" frame of the bound state, defined by 
P = Pi+p 2 = 0) that the square M 2 of any bound-state mass M is real, respectively. 

• In the center-of-momentum frame of the bound state, any solution $(p) of Salpeter's 
equation (1) is subject to the relation (cf. Eq. (2.17) of Ref. [8] or Eq. (11) of Ref. [10]) 



M||$| 



-/ 

r &'p r & J t 

+ J ThxYJ (27 



Tr 



+ 



^[ El {p)+E 2 {p)]Ti[^{p)^p 

d 3 p r d 3 q 
(2^)3 J (2^p 
d 3 p 
(2^)3 

d 3 p f d 3 q 
J^f! (27p 



$ f (p) 70 [K(p, q) &(q)\ 7o 
[£\(p)+£ 2 (p)]Tr[$t(p)$(p) 

]T V T (p, q) Tr [&(p) 7o T $( 9 ) T 7o ] . (15) 



For all interaction kernels such that / d 3 p / d 3 q Tr &(p) 70 [K{p, q) $(q)] 70 is real, 



d 3 p f d 3 q 



(27r)3 y i27r)3 {^[^(p)lo[K(p,q)Hq)ho]) 
d 3 p f d 3 q 



7o 



the right-hand side of Eq. (15) is, of course, real. In this case, remembering the reality 
of the square of the Salpeter "norm" , the mass eigenvalues M of all bound states with 
nonvanishing "norm" ||$|| of the associated Salpeter amplitudes $, that is, ||$|| 2 7^ 0, 
are real: M* = M [8,23]. The needed reality of the interaction term in Eq. (15) holds, 
for instance, for all Bethe-Salpeter kernels which are sums of terms of Dirac structure 
r®r such that the Dirac matrices T = 7o T are (anti-) Hermitian, = ±f, with the 
associated Lorentz-scalar potential function Vr(p, q) satisfying V^(q, p) = Vr(p, q). 1 



1 Thus, (at least) for the present class of interactions any solution of Salpeter's equation (1) with nonreal 
and therefore purely imaginary mass eigenvalue M* 7^ M = ±i \M\ must have vanishing "norm" || $|| 2 = 0. 
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• The comparison of Salpeter's equation (1) with its Hermitian conjugate (sandwiched 
between 7 from the left and 70 from the right) reveals [8,23] that the solutions of this 
eigenvalue equation always occur in pairs of the form ($, M) and (\l/, — M*), with the 
Salpeter amplitude ^(p) related to the Hermitian conjugate of its "partner solution" 
provided the interaction kernel satisfies [K(p, q) ^{q)}^ = 70 [K(p, q) ^(q)] 70; 
the squares of the "norms" of conjugate solutions have opposite sign: \\^\\ 2 = — \\&\\ 2 - 
For the case of equal masses m\ = m 2 of the bound-state constituents [8] , the relation 
between <3> and \l/ is particularly simple, viz., = 70 &(p) 7o- For the general case 
of unequal masses mi 7^ m 2 of the bound-state constituents [23] the relation of $ and 
^ involves, in addition, the charge-conjugation matrix C defined, e.g., by C = i 7 2 7 . 
This doubling of the eigenvalues to pairs of opposite sign is well known from solutions 
to the RPA equation [23,37,38]. The requirement on the interaction kernel holds, for 
example, for every Bethe-Salpeter kernel which is a sum of terms each of which is the 
product of a real Lorentz-scalar interaction function Vp(p, q) = V^(p, q) and a Dirac 
structure T®T such that all Dirac matrices T = 70 T are (anti-) Hermitian, Y* = ±r. 
Furthermore, the relation between $(p) and holding within the particular class 
of interactions defined by the requirement [K(p, q) = 70 [K{p, q) ^(q)} 70 also 

entails that a Salpeter amplitude $ solving the Salpeter equation for a nondegenerate 
zero mass eigenvalue M = necessarily has vanishing Salpeter "norm:" ||$|| =0 [23]. 

In particular, any momentum-space interaction function Vr(p, q) which can be represented 
as the Fourier transform Vr(p, q) = Vr(p—q) = J d 3 pexp[— i (p— q)-x] Vp(r) of some (real) 
central potential Vr(r) = V^(r) in configuration space experienced by the two bound-state 
constituents at distance r = |cc| satisfies both V^(p, q) = Vr(p, q) and Vr(q,p) = Vr(p, q). 
Thus, for Dirac matrices V in the kernel K(p, q) = J2r ^r(p, q) obeying 70 70 = ±r 
both conditions on the kernel K(p, q) required for the above-mentioned spectral properties 
to hold are simultaneously fulfilled. In this case, the spectrum of mass eigenvalues M of the 
Salpeter equation (1) can only involve, in the complex-M plane, pairs (M, — M) of opposite 
sign on the real axis (excluding the origin), with the associated Salpeter amplitudes having 
Salpeter "norms" squared of opposite sign, and/or points M = —M* on the imaginary axis 
(including the origin), with corresponding Salpeter solutions of vanishing Salpeter "norm." 

Merely for illustrative purposes let us now generalize, for a few examples, the treatment 
applied to the reduced Salpeter equation (4) in Sec. 3 through Sec. 6 to the case of the (full) 
Salpeter equation (1). As recalled in Sec. 3, Salpeter amplitudes $(p) for bound states with 
spectroscopic label n 1 Jj and, therefore, all Salpeter amplitudes <3>(p) for bound states with 
spin-parity-charge conjugation assignment J PC = CT + consist of two independent Salpeter 
components <pi(p) and 02 (p)- The corresponding Salpeter equation will therefore reduce to 
a system of two coupled (integral) equations for two radial wave functions <pi(p) and 02 (p)- 2 

In order to analyze simultaneously Bethe-Salpeter interaction kernels of Lorentz-scalar 
and time-component Lorentz- vector Dirac structure, we introduce a simple sign factor o by 




+ 1 for T ® T = 7 ® 7 (time-component Lorentz- vector interactions) , 
— 1 for r <g> T = 1 ® 1 (Lorentz-scalar interactions) . 



2 As recalled in Sec. 3, for pseudoscalar bound states the constraints on the solutions &(p) of the reduced 
Salpeter equation (4) imply <pi (p) — <p 2 (p) = <p{p)- Consequently, for a given Lorentz structure r® T of the 
interaction kernel the radial eigenvalue equation representing the reduced Salpeter equation can be derived 
by adding the two radial eigenvalue equations related to Salpeter's equation (1) after letting <j)\{p) = (f>2{p)- 
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(Lorentz-scalar and time-component Lorentz- vector confining interaction kernels and their 
linear combinations attracted, for purely phenomenological reasons, particular attention in 
the Bethe-Salpeter descriptions [8-15,39-41] of hadrons as bound states of quarks [42,43].) 
With the help of the parameter a and the interaction functions Vl(p, q) (L = 0, 1) of Sec. 4, 
the radial remnants of the Salpeter equation (1) for pseudoscalar bound states are, for both 
Lorentz-scalar (r®r = 1(8)1) and time-component Lorentz-vector (r®r = 7°(87°) kernels, 



2E(p) 2 (p) + J 



dq q 2 

WT 2 



aV (p, q )Mq) = MMp) , 



rl 2 



m 



E(p)E(q) 



V (p,q) + 



pq 



E(p)E(q) 



Vi(p,q) 



Mq) = M 2 (p) 



This set of equations exhibits a very peculiar structure: For M — 0, the equations decouple. 
For M^O, from one of the equations one of the two independent Salpeter components, say 
<j>i(p), may be expressed in terms of the other, 2 (p), and inserted into the other equation in 
order to reformulate this eigenvalue problem for the square M 2 of the bound-state mass M: 

dq q 



4E 2 (p)Mp) + 2E(p) J ^LaV (p,q)M<l) 



+ 



/ 



dqq 2 



E(p)J (2n) 2 



crm 2 V Q (p,q)+pqV 1 (p,q) (j) 2 (q) 



+ 



dq q 2 

o v ' 
M 2 2 (p) 



m 



a 



E(p)E(q) 



V (p, q) + 



pq 



E(p) E(q) 



Vi(p,q) 



dkk 2 
(2^ 



aV (q,k)(j)2(k) 



Following Sec. 5, specifying the interactions to harmonic-oscillator form allows to trade any 
integration over the potential functions V L (p, g), L — 0, 1, for the differential operator 



4 E\p) -^- (a(p 2 + 2 m 2 ) + p £><°> p - 2) 



E(p) 



+ 



W)^ Df)+apDf)p -^m Df) 



i 2 ( p ) = m 2 Up) ■ 



This is still a homogeneous linear differential equation but of fourth order. Thus its analysis 
requires techniques or tools beyond those familiar from the study of Schrodinger operators. 

The Salpeter equation for a Bethe-Salpeter kernel of the BJK [16,17] Lorentz structure 
T(E)T = | (7^ (8 7 M + 75 (8 75 — 1 (8 1) constitutes an exceptional case; there all interactions 
enter only in one of the two relations which form the set of coupled radial integral equations 
describing pseudoscalar bound states, the other relation being of "merely" algebraic nature: 



2 + 2 

2E( p )Mp) = mMp) ■ 



V (p,q) M<1) = MMP) 
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Merging these two relations generates an eigenvalue problem for M 2 , posed equivalently by 

" J 2 

aq q 



AE\p) 0i(p)+4 



(27T) 



2 H(p,g)^(g)0i(g) = M 2 0x(p) , 



AE 2 (p) <f> 2 {p)+±E{p) 



dqq 2 



V (p,q) <h(Q) 



M 2 . 



ar 



For the harmonic-oscillator interaction V(r) 
virtue of Eq. (5), i.e., q 2 Vo(p, q) = —(2tt) 2 a 5(p—q), the ordinary differential equations 



2 these eigenvalue equations become, by 



E 2 (p)-aD^E(p)]<p 1 (p) = M 2 MP), 
E 2 (p)-aE(p)DfAUv) = M 2 MP)- 



Applying the substitution E{p) 4>\{p) oc 02 (p) proves the equivalence of these formulations. 
In contrast to the case of an arbitrary interaction kernel, for the B JK Lorentz structure our 
harmonic-oscillator Salpeter problem reduces to a single second-order differential equation. 
The differential operator on the left-hand side of this eigenvalue problem is not self-adjoint. 
Nevertheless, according to the general properties of eigenvalues of the Salpeter equation (1) 
summarized at the beginning of this section, the spectrum of eigenvalue squares M 2 is real. 
Recalling, for a > 0, our reasoning of Subsec. 8.1 for H v = — A+U with auxiliary potential 

E(p) M 2 



U(p) = 



a AaE(p) ' 



it is trivial to demonstrate, by similar arguments, that the spectrum of squared eigenvalues 
M 2 and, as a consequence thereof, the spectrum of mass eigenvalues M are purely discrete. 



10 Summary, Conclusions, and Outlook 

The present investigation has been devoted to an exploration of the conditions under which 
the reduced Salpeter equation with confining interactions has stable bound-state solutions. 
For harmonic-oscillator interactions, the reduced Salpeter equation becomes in momentum 
space either an algebraic relation or a second-order ordinary differential equation involving 
the Laplacian A acting on states of angular momentum i = (i.e., our differential operator 
introduced in Sec. 5). For real harmonic-oscillator couplings, all corresponding spectra 
are real. For pseudoscalar states, where instabilities are expected to appear first, we showed 
that, depending on the Lorentz nature of the kernel, the resulting spectrum is either purely 
continuous or entirely discrete, consisting of isolated mass eigenvalues of finite multiplicity. 

As a by-product, the same analysis proves the boundedness from below of all the spectra 
for appropriate choice of the respective sign of the harmonic-oscillator coupling constant a. 

• For interaction kernels of time-component Lorentz- vector structure (r®r = 7°®7°), 
Lorentz-pseudoscalar structure (r®r = 75(^75), the (eventually simple) BJK [16,17] 
structure, T®r = | (7^ ®7 M +75(^75 — 1<8>1), and, if m — 0, of Lorentz-scalar structure 
(r®r = 1(8)1) or Lorentz-vector structure (rtgT = 7^(87^), each harmonic-oscillator 
ordinary differential equation of Sec. 6 is, for arbitrary a > 0, the eigenvalue equation 
of some positive operator. The entire spectrum of any such operator must be positive. 
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• For an interaction kernel of Lorentz-pseudoscalar structure = 7 5 Cg>7 5 ), and any 
m^fl, Eq. (9) is, for finite a < 0, the eigenvalue equation of an operator which is not 
positive but bounded from below, by 2m+3a/(2m 2 ), as is its (continuous) spectrum. 

• For m^O, the ordinary differential equations, in Sec. 6, related to interaction kernels 
of Lorentz-scalar structure, r<g>r = l(g)l, or Lorentz-vector structure, r<g>r = 7 M £g>7 M , 
are not of the (standard) form of eigenvalue equations of some Schrodinger operators. 
In order to decide on the nature of their spectra we have to rely on the transformation 
of Subsec. 6.1 to Schrodinger-like auxiliary operators. For a < in the Lorentz-scalar 
case and a > in the Lorentz-vector case, all eigenvalues £»(M), i — 0, 1, 2, . . . , of the 
latter operators can be shown to be entirely discrete for all M, and strictly decreasing 
functions of M. Because of their strict decrease, with increasing bound-state mass M, 
the zeros of the lowest trajectories Sq(M) define not necessarily positive lower bounds 
on the spectra of mass eigenvalues M, proving both spectra to be bounded from below. 

Altogether this provides a rigorous proof of the stability of the considered bound states: 
their energies form (for couplings of suitable sign) real discrete spectra bounded from below. 

Our findings point in the same direction as the observations made in a purely numerical 
analysis [44] aiming at the description of quark-antiquark bound states with the help of the 
full Salpeter equation (1). In this investigation the (instantaneous) interaction between the 
bound-state constituents is modeled by the sum of a Coulomb-type short-range interaction 
(arising from one-gluon exchange between quark and antiquark) and a rather sophisticated 
confining interaction, which interpolates between a harmonic-oscillator-type behaviour for 
small inter-quark distances r, or small masses of the bound-state constituents, and a linear 
rise for large inter-quark distances r, or large masses of the bound-state constituents. From 
a point of view of principle it is a pity that this study has been performed merely at a single 
point in free-parameter space, determined (for an equal-weight mixture of time-component 
Lorentz-vector and Lorentz-scalar Dirac structures of this particular confining interaction) 
from a fit of the mass spectrum of experimentally observed mesons. The authors of Ref. [44] 
arrive at the following conclusions (for bound states composed of equal-mass constituents): 

• For bound states of heavy-mass constituents, the solutions prove to be stable for both 
time-component Lorentz-vector and Lorentz-scalar Dirac structures of the confining 
interaction and, consequently, for any linear combination of these Lorentz structures. 

• For bound states of light-mass constituents, the solutions are still all stable for a pure 
time-component Lorentz-vector Dirac structure of the confining interaction but turn 
out to be mostly unstable for a pure Lorentz-scalar spin structure of the confinement. 

However, in Ref. [44] it was also (numerically) shown that neglecting the "negative-energy" 
components of the Salpeter amplitude $(p) by simply imposing onto $(p) the requirement 
Aj"(p 1 ) $(p) A2"(p 2 ) = (almost totally) removes the instability of the bound-state energy 
levels otherwise showing up in the case of a Lorentz-scalar confining interaction. According 
to our discussion of the component contents of all solutions of both the full and the reduced 
Salpeter equations with respect to the energy projectors Af (p) in Sec. 2, imposition of such 
constraint on $(p) is tantamount to the consideration of the reduced Salpeter equation (4). 

Our approach can be clearly generalized [35] to analyze not only Salpeter 's equation (1) 
but, in a similar way, three-dimensional reductions of the Bethe-Salpeter equation different 
from the Salpeter equation; some of these reductions are reviewed, for instance, in Ref. [45]. 
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The intention of the present study was to perform, as a purely theoretical investigation, 
a rigorous analysis of the spectral properties of the Salpeter equation (4), irrespective of its 
four-dimensional origin within the Bethe-Salpeter formalism. Our pragmatic point of view 
seems to be justified by the fact that in the past exactly this bound-state equation has been 
applied in numerous treatments of, for instance, quark-antiquark bound states in quantum 
chromodynamics (for details, consult the reviews in Refs. [42,43]). As already mentioned in 
the Introduction, our interest in the problem of the particular type of instabilities discussed 
here has been aroused by their occasional observation [10,44], and the attempts to arrive at 
full analytic understanding of them [39-41]. Nevertheless, we feel obliged to add a few brief 
comments on the practical relevance of this framework and the implications of our findings. 

1. From the very beginning, our considerations have been confined to three-dimensional 
reductions of the Bethe-Salpeter formalism, obtained by assuming all interactions to 
be instantaneous in the center-of-momentum frame of the two-particle systems under 
study. The relevance of all statements about the presence or absence of instabilities in 
any three-dimensional reduction (as given by Salpeter's equation) for the situation in 
the four- dimensional Bethe-Salpeter formalism remains unclear; its discussion would 
necessitate a spectral analysis of its own, which will become much more involved than 
our analysis reported above, but is well beyond the scope of the present investigation. 

2. As recalled in Sec. 2, in addition to some instantaneous approximation any derivation 
of, in particular, the Salpeter equation relies on the assumption of free propagation of 
all bound-state constituents, with a constant (i.e., momentum-independent) effective 
or constituent mass which should adequately parametrize some significant part of the 
dynamical self-energy effects. Needless to say, the free-propagator assumption cannot 
be compatible with a confining interaction and is therefore conceptually problematic: 
In quantum field theory, the Dyson-Schwinger equations relate the propagators, that 
is, the two-point Green functions, to the n-point Green functions which represent the 
interactions in the Bethe-Salpeter equation. Therefore, propagators and interactions 
cannot be chosen independently from each other. Thus in non-Abelian gauge theories 
the simultaneous assumption of free propagation of the bound particles and confining 
interactions (as induced by quantum chromodynamics) are intrinsically inconsistent. 
Moreover, phenomena such as the dynamical breakdown of chiral symmetry can only 
be taken into account by retaining the exact propagators, obtained as the solutions of 
the Dyson-Schwinger equation for the corresponding two-point Green function of the 
bound-state constituents. Their proper incorporation is crucial for the interpretation 
of the lowest-lying pseudoscalar quark-antiquark bound states: as Goldstone bosons. 
To make a long story short, it is, of course, very desirable to have also in one's favorite 
instantaneous bound-state equation the exact fermion propagators at one's disposal. 

One of the attempts to retain in some three-dimensional reduction (as far as possible) 
the original wave-function renormalization and mass functions parametrizing the full 
propagators of both fermionic bound-state constituents resulted in the instantaneous 
bound-state equation proposed in Ref . [46] . The first tentative exploration of some of 
the implications of this generalization [47] of Salpeter's equation for quark-antiquark 
bound states may be found in Ref. [48] . The latter study utilizes exact propagators of 
light quarks, as extracted within a "renormalization-group-improved rainbow-ladder 
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truncation" (which scheme has the undeniable advantage to preserve the axial- vector 
Ward-Takahashi identity) applied to both the quark Dyson-Schwinger equation and 
the meson Bethe-Salpeter equation [49]. Given the formalism introduced in Ref. [46], 
our (obvious) next step is a similar study [34,50] for our full-propagator version of the 
Salpeter equation (4). Our preliminary results indicate that for reasonable behaviour 
of both bound-state constituents' propagator functions, i.e., nontrivial wave-function 
renormalization and dynamical mass, stability (in our sense) can be achieved [34,50]. 
The significance of these findings for the four-dimensional Bethe-Salpeter formalism, 
which includes all the constituents' self-energy effects, may be judged by future work. 

3. Of course, one might argue that a sufficiently precise purely numerical solution of the 
bound-state equation in use may suffice to settle the stability issues once and forever. 
However, as demonstrated, e.g., by the not really conclusive findings of Refs. [10,44] a 
numerical analysis may give but a hint of potential problems with unstable solutions. 
Thus, it is our conviction that a genuine understanding of the origin of these troubles 
and a compelling solution of this problem may be gained only by some analytic proof. 

Hence, as our conclusion let us stress once again that, for the solutions of the instantaneous 
Bethe-Salpeter formalism, at least, the kind of analysis proposed above provides a rigorous 
answer to the question of stability, defined by our requirement that each bound state found 
belongs to a real and discrete spectrum that is bounded from below. Note that the requested 
reality of the bound-state masses precludes, for example, any solutions of tachyonic nature. 
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